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We estimate the probability of detecting gamma-rays from the annihilation of neutralino dark 
matter in dense, central regions of Milky Way substructure. We characterize Galactic substructure 
statistically based on Monte Carlo realizations of the formation of a Milky Way-like halo using a semi- 
analytic method that has been calibrated against the results of high-resolution N-body simulations. 
We find that it may be possible for the upcoming experiments GLAST and VERITAS, working in 
concert, to detect gamma-rays from dark matter substructure if the neutralino is relatively light 
(M-^^ 100 GeV), while for 500 GeV such a detection would be unlikely. We perform most 

of our calculations within the framework of the standard ACDM cosmological model; however, we 
also investigate the robustness of our results to various assumptions and find that the probability of 
detection is sensitive to poorly-constrained input parameters, particularly those that characterize the 
primordial power spectrum. Specifically, the best-fitting power spectrum of the WMAP team, with 
a running spectral index, predicts roughly a factor of fifty fewer detectable subhalos compared to the 
standard ACDM cosmological model with scale-invariant power spectrum. We conclude that the lack 
of a detected gamma-ray signal gives very little information about the supersymmetric parameter 
space due to uncertainties associated with both the properties of substructure and cosmological 
parameters. 

PACS numbers: 95.35.-|-d, 98.62.Gq, 98.35.Gi, 14.80.Ly 



I. INTRODUCTION 



A standard cosmology (ACDM) has emerged in which 
the Universe is spatially flat and its energy budget is bal- 
anced by ~ 4% baryonic matter, ~ 26% cold, collision- 
less dark matter (CDM), and roughly 70% dark energy 
or a cosmological constant (A) Q, |^. The growth 
of structure is seeded by density fluctuations supposedly 
generated during an early epoch of inflation 4] . The pri- 
mordial power spectrum of fluctuations is expected to 
have a nearly scale-invariant form, P{k) oc /c", n ~ 1. 
CDM dominates the baryons in the matter budget, and 
luminous galaxies form within halos of CDM, where the 
dark matter potential wells trap and compress baryons 
so that they may cool and condense. In this paradigm, 
structure forms hierarchically: small mass objects col- 
lapse first and merge into larger objects over time f^. 
Many of the small halos subsumed by objects that grow 
to become present day galactic halos survive as gravita- 
tionally self-bound substructure or "subhalos." In fact, 
CDM theory predicts that halos similar to that of the 
Milky Way (MW) should host hundreds of subhalos that 
are apparently about the same size as the observed satel- 
lite galaxies of the MW 6]; however, there are only eleven 
known satellite galaxies within ~ 300 kpc of the MW 
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(seeRef. for a census of the Local Group) . This dis- 
crepancy has been dubbed the "dwarf satellite problem" 
(DSP). 

The robustness of the DSP is uncertain and some au- 
thors have argued that the observed MW satellites may 
correspond directly to the eleven or so most massive sub- 
halos observed in simulations Alternatively, others 
have proposed solutions to the DSP that fall into two 
very broad categories. In the first, the number of sub- 
halos of the appropriate size and density structure is re- 
duced such that all such subhalos host luminous satellite 
galaxies. This outcome can be achieved either by chang- 
ing the properties of the dark matter Q or by modifying 
the spectrum of density fluctuations that seed structure 
growth In the second class of solutions, MW- 

like halos do play host to hundreds of subhalos of the size 
expected to host its dwarf satellites, but luminous galax- 
ies do not form in ~ 90% of these subhalos and they 
remain dark and undetectable. In this case, baryons do 
not cool and condense in subhalos, perhaps because they 
are heated by the photo-ionizing background ,13j|] , by 
supernova feedback ll4|, or some other feedback mech- 
anism. Efforts to detect non-luminous substructure in 
galactic halos through lensing effects 15], tidal streams 
|l6j| , or any other method therefore promise to distinguish 
between these alternatives. As such, these observations 
may provide a crucial test of the ACDM paradigm and 
reveal important information about structure growth and 
galaxy formation. 

If the standard CDM paradigm is correct, then galac- 
tic halos should be teeming with substructure and if the 
dark matter is in the form of a supersymmetric particle 
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species produced in the early Universe, then it is possi- 
ble for this substructure to be lit up by the annihilations 
of the dark matter particles into gamma-rays. Super- 
symmetry (SUSY) is the most popular extension of the 
Standard Model (SM) of particle physics (for a review of 
SUSY, see Ref. with several intriguing properties, 

one of the most cosmologically significant of which is that 
it admits natural CDM candidates. In the most popular 
models, especially from the standpoint of building SUSY 
grand unified theories, the conservation of R-parity guar- 
antees that the lightest supersymmetric particle (LSP) is 
stable. Moreover, there is a large swath of SUSY param- 
eter space in which the LSP is weakly-interacting and 
may be produced in the early Universe with an appro- 
priate relic density to serve as the CDM. In the simplest 
models that meet accelerator constraints, the LSP is the 
lightest neutralino (x), or the lightest mass eigenstate 
formed from the superposition of the two CP-even Hig- 
gsinos, the W'^ino, and the Bino (for a review of SUSY 
that focuses on SUSY dark matter, see Ref. 

Several authors have considered the possibility of di- 
rect detection of neutralino dark matter from the gamma- 
ray flux produced by neutralino annihilations in the 
Galactic center (T^. Assuming that the Milky Way re- 
sides within a standard NFW-like halo, the flux from 
the Galactic center may be considerably larger than that 
from substructure (by perhaps as much as a factor of 
~ 10'^) 125; however, there are large uncertainties in 
the dark matter density structure of Galaxy-sized ha- 
los, largely due to the unknown effect of baryons on 
the dark matter distribution [23|. In order to mitigate 
these uncertainties and to test the CDM paradigm, it 
may be advantageous to consider annihilations in CDM 
substructure. Other authors have considered this pos- 
sibility H m m m m El m, usuaUy with an eye 
toward constraining the minimal supersymmetric stan- 
dard model (MSSM) parameter space. Generally, these 
studies have used overly optimistic prescriptions for the 
spatial distribution and density profiles of CDM subha- 
los that are not supported by detailed studies of structure 
formation in ACDM cosmologies. 

In this paper, we study neutralino annihilation in halo 
substructure using a rather different approach. We use 
a semi-analytic method to estimate the properties of the 
subhalo population. Past studies of neutralino annihila- 
tion in substructure n m m m m did not account 
for the fact that the density structure of the subhalos and 
the spatial distribution of substructure of a given mass 
depend upon their accretion histories. The semi-analytic 
model we use allows us to account for these correlations, 
albeit in an approximate way. In addition, the model 
allows us to generate statistically significant results for 
a variety of input paramters by examining a large num- 
ber of realizations of MW-like halos. Thus we can esti- 
mate the likelihood of observing gamma-rays from neu- 
tralino annihilations in dark subhalos. In this sense, our 
approach is aimed more toward serving as a guide for 
observations and is complementary to the recent numer- 



ical study of Stoehr et al. We demonstrate that 

observations of this kind will likely not yield meaningful 
constraints on SUSY and explicitly show how predictions 
for the expected flux and the number of observable sub- 
halos are rather sensitive to uncertainties regarding the 
background cosmology. In particular, we show that the 
likelihood of directly observing subhalos via neutralino 
annihilations is very sensitive to the initial power spec- 
trum of density perturbations on sub-galactic scales. 

The outline of this manuscript is as follows. In Sec- 
tion m we review the properties of dark matter halos 
observed in N-body simulations and discuss our model- 
ing of halo substructure. In Section UTTI we calculate the 
expected flux from CDM subhalos, describe the various 
backgrounds that this signal must overcome and outline 
conditions for detectability. In Section llVl we present our 
results on the detectability of substructure via gamma- 
rays from neutralino annihilations. Lastly, we summa- 
rize and draw conclusions from this work in Section |V| 
Throughout most of this work, we assume the standard 
ACDM cosmological model (e.g., [IQ) with Qm = 0.3, 
r^A = 1 - l^M = 0.7, 17b/i^ = 0.02, h = 0.72, and 
a standard, scale-invariant primordial power spectrum 
of initial density fluctuations, -P(fc) oc fc" with n = 1. 
We also consider the consequences of adopting a primor- 
dial power spectrum with a running power law index 
dn/d\nk = —0.03, normalized to as — 0.84 as advo- 
cated by the recent analysis of the Wilkinson Microwave 
Anisotropy Probe (WMAP) team 0. 



II. THE STRUCTURE AND SUBSTRUCTURE 
OF DARK MATTER HALOS 

In this Section, we discuss the density structure of 
dark matter halos as well as the properties and charac- 
teristics of halo substructure based on our semi-analytic 
model. The size of a halo can be quantified by its 
virial Mass Mv, virial radius Rv, or equivalently its virial 
velocity = GAU/Rv The virial radius is defined 
as the radius within which the mean density is Avir 
times the mean matter density of the universe pM, so 
that Mv = 4:7rpM^vii{z)Rl/'i- The virial overdensity 
can be estimated using the approximation of spherical 
top-hat collapse. We compute Avir (2) using the ap- 
proximate fitting formula of Bryan and Norman |28j : 
Avh.(2) ^ (IStt^ -t- 82a; - 39a;2)/f7M(z), where x + 1 = 
nuiz) = nuil + z)^/[nM{l + z)^ + f^A] and nuiz) is 
the ratio of the mean matter density to critical density 
at redshift z. In the ACDM cosmology that we adopt, 
Avir(z = 0) ~ 337 and Avir{z) 178 at high redshift, 
approaching the standard CDM {Qm — 1) value. For ref- 
erence, the virial radius can be written in terms of Mv 
asi?v = 13.6/i-ikpcAfg/^[f7MAvir(z)/337]-i/3(l + z)^i 
where Mg = Mv/lO^ h'^ Mq. 

The matter density profiles of CDM halos have been 
studied extensively in numerical simulations. The 
spherically-averaged density profile proposed by Navarro, 
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Frenk, and White ^ (hereafter NFW) 

,W=..(^)"(i + ^)". (1) 

seems to represent best the structure of halos in the most 
recent studies of high-resolution numerical simulations 
|30| and seems to be further buttressed by other theoret- 
ical considerations [sil Is^ . These studies rule out the 
singular isothermal sphere with p{r) oc and begin 
to disfavor the very steep profile with p(r) cx r~^'^ pro- 
posed by Moore et al. [33 as faithful representations of 
the predictions of CDM on small scales. In fact, Power 
et al. jlOl] observe that the inner profiles of halos be- 
come increasingly shallow with decreasing radius all the 
way down to the minimum radius at which their conver- 
gence criteria are met. Therefore, even the assumption 
of the NFW profile is likely to overestimate gamma-ray 
flux from dark matter annihilation because this flux is 
very sensitive to the dark matter density at radii below 
the resolution limits of present N-body simulations. A 
definitive resolution of this issue awaits further numer- 
ical work. In most of this study, we assume the NFW 
profile with inner power law p cx r^^, but we also explore 
the effect of adopting the steeper profile of Moore et al. 
[33j | on our results. 

The relative concentration of an NFW halo of a given 
mass is determined by the NFW scale radius Tg or equiv- 
alently, the NFW concentration parameter Cv = Rv/^a- 
This is related to the characteristic NFW density by 
Ps = pMAvirC^/3/(cv), where /(x) = ln(H-a;)-a::/(H-a;). 
Numerical studies have shown that the concentrations of 
dark matter halos are set by their mass accretion histories 
UMI^HEHS' In particular, the concentration param- 
eter of a halo of a given mass is strongly correlated with 
a suitably-defined epoch of formation for typical halos 
of that mass (see Refs. [2^ |3^ for details). The 
important point is that the structural parameters of the 
NFW profiles that describe halo substructure are essen- 
tially fixed by the choice of cosmology and are not set 
by, for example, requiring a subhalo of a given mass at 
a given galacto-centric radius to be compact enough to 
resist tidal stripping as has been supposed in previous 
studies (e.g., as in Ref. jl^)- Moreover, because of this 
correlation between subhalo concentration and halo for- 
mation epoch, the structural parameters of halos depend 
strongly upon cosmological parameters such as the lin- 
ear, rms amplitude of density fluctuations smoothed over 
spheres of radius 8 h~^Mpc, as, and the power law in- 
dex of the primordial power spectrum, n 37] (of course, if 
the normalization of the power spectrum is fixed on large 
scales by measurements of the cosmic microwave back- 
ground anisotropy, for instance, a change in n implies a 
specific shift in erg). As we discuss below, the gamma- 
ray flux from CDM annihilations in substructure is quite 
sensitive to the density of the inner regions of halos. 
This implies that the probability of detecting gamma- 
ray flux from neutralino annihilations in substructure can 
be a strong function of cosmological parameters, even in 



highly-idealized calculations like those that we present 
here. 

In order to account for the correlations between halo 
mass, redshift, and concentration we adopt the simple, 
semi-analytic prescription of Bullock et al. ^5|| (BOl) in 
our modeling of substructure. We use this model to set 
the concentrations of all substructure halos as described 
below. The BOl model reproduces the mean Cv — 
relation observed in N-body simulations and provides an 
estimate of the statistical scatter. The model has been 
tested successfully against standard ACDM simulations, 
tilted ACDM simulations, simulations of so-called "stan- 
dard" CDM (i.e., = 1), scale- free power law CDM 
models ,3^ ,33 , and warm dark matter cosmologies [s^ 
for halos in the mass range 10^^ Mv/Mq^ 10^'*. In ad- 
dition, recent results from new simulation data together 
with a re-analysis of older simulations suggest that the 
BOl model accurately represents the mean Cv — rela- 
tion all the way down to ~ 10^ M^p, and also supports 
the universality of the NFW profile |4a.l4ll |. 

The first step in our calculation is to estimate the spa- 
tial distribution and structural properties of substruc- 
ture in the MW halo. The properties of substructure in 
CDM halos are determined through an endless compe- 
tition between subhalo accretion and destruction due to 
tidal forces and dynamical friction. Several authors have 
studied substructure distributions in simulated CDM ha- 
los 6, 42, 43, 44], typically in cluster-sized halos. In 
this study, we adopt a complementary, semi-analytic ap- 
proach to estimate the properties of the substructure 
population as described in Ref. (hereafter ZB03). 

We give a brief summary of the model below, see Ref. 
|Tll | for further detail. 

First, we generate the merger histories of MW-sized 
host halos using the extended Press-Schechter formalism 
[isf . In particular, we employ a modified version of the 
merger tree algorithm proposed by Somerville and Kolatt 
|46| . This allows us to generate a list of masses and accre- 
tion redshifts for all halos above a given mass threshold 
-^^min, that merged into the present day MW-like halo 
during its formation. Second, we assign each subhalo an 
NFW concentration at the accretion event. To do this, 
we use the Bullock et al. 35] model to calculate the 
mean value of the logarithm of the concentration param- 
eter, (log(cv)), for the mass of each subhalo at the time 
of accretion. We then compute the actual value of Cv 
that we assign to each subhalo by selecting Cv randomly 
from a log-normal distribution with standard deviation 
(T(log(cv)) = 0.14. Bullock et al. [3^ determined this 
to be a good approximation for the statistical scatter of 
Cv at fixed mass in their simulations. Finally, we track 
the subhalo's orbit in the potential of the host from the 
time of accretion until today, in a manner similar to that 
of Taylor and Babul j47l] , in order to determine whether 
or not the subhalo is destroyed by tidal forces once in- 
corporated into the parent halo. Using this model we 
construct 100 statistical realizations of a MW-sized halo 
with Afv = 1.4x 10^^ M0, using a minimum subhalo mass 
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of Mmin = 10^ Mq, and ten realizations of a MW-sized 
halo with M^in = 10^ Mq. 

This technique yields an estimate of the substructure 
population of MW-sized halos that is in approximate 
agreement with the radial distributions, mass functions, 
and velocity functions of high-resolution N-body simula- 
tions pdj . This method offers improvement over some 
previous estimates of substructure populations in stud- 
ies of neutralino annihilation in substructure because we 
account for the known correlations between mass accre- 
tion history and the structure of both the host and its 
subhalos and we are able to study a statistically signif- 
icant number of hosts. Like previous numerical studies 
E3j ) this method reveals a core-like behavior in 
the radial distribution of subhalos themselves for radii 
30 — 70 kpc, with the exact value depending upon the 
subhalo mass cut one considers. This is the well-known 
result that subhalos are anti-biased with respect to the 
dark matter due to tidal effects. The distribution of sub- 
halos only follows the p oc r^"^ behavior of the mean dark 
matter distribution at fairly large halo-centric radii. In 
addition, our method complements N-body work. First, 
it suffers from no inherent resolution effects such as the 
issue of "overmerging," which may cause simulations to 
underestimate the amount of substructure near the cen- 
ters of Galaxy-sized halos 42]. If resolution issues are 
not a problem in numerical simulations, then this method 
provides an overestimate of gamma-ray flux from subha- 
los. This is a conservative approach in the context of 
this paper. Second, our method allows us to generate a 
large number of MW halo realizations so that we may 
estimate the probability of observing gamma-rays from 
neutralino annihilations in otherwise dark substructure, 
given the scatter in the subhalo distribution in any given 
host, and demonstrate the effects of cosmological param- 
eters on the probability of detection. 

With the gross properties of the subhalo distribution in 
place, another issue that must be dealt with is the distri- 
bution of matter in the very central regions of subhalos. 
This is where most of the luminosity due to annihila- 
tions originates because the annihilation rate scales as 
the square of the mass density (see below). There will 
be some limit to the central density of the halo. The 
densities achieved in the inner regions of halos may be 
limited simply due to the fact that neutralinos in very 
dense regions will annihilate, so something closer to a 
nearly constant density core might be expected. In this 
case, the core size is set by a competition between mass 
in-fall and neutralino annihilation, so a simple approxi- 
mation for the core radius r^o follows from equating the 
annihilation rate to the rate of in-falling material, 

MJ p,{r,,o){cj\v\) ^ ^Rl/GM,, (2) 

and solving for rc.o- In Eq. ([SJl, is the mass of 
the neutralino and (<t|w|) is the thermally-averaged total 
cross section times relative velocity for neutralino annihi- 
lation. We assume that the density at all points interior 
to Tcfi is given by the NFW profile density at that point. 
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FIG. 1: Left: The relationship between plri (this quan- 
tity is proportional to subhalo luminosity, see Section lili Al 
and Section ITVIl and the mass M, of Galactic subhalos. The 
points represent a scatter plot of these quantities for all sur- 
viving halos in one realization of our substructure model. The 
thick solid line represents the values of ^ a function 

of mass using the prescription of Ref. 1221 for the proper- 
ties of the subhalos that they place nearest the sun. The 
dashed lines represent the 99.9 percentile range of subhalo 
properties in our model. Right: The relationship between 
the quantity pi rj and the distance between the sun and the 
subhalo, d. The points represent a scatter plot of pi r'l vs. 
d for all of the surviving subhalos with masses in the range 
10^ < M/Mq < 5 X 10^ in ten model reahzations. The thick 
solid line represents the range of the pi rl-d relationship for 
subhalos in the same mass range obtained by setting the prop- 
erties of the nearest subhalos according to the prescription of 
[2^ . Notice that our modeling suggests that subhalos have 
significantly lower luminosity (proportional to plrl) and are 
typically further away than what was assumed in Ref. |2^ . 



Pc = p{fc,o)- In actuality, most of the gamma-ray lumi- 
nosity arises from the inner core regions, which are orders 
of magnitude smaller than has been reliably probed by 
N-body simulations, so there is little justification for as- 
sumin g th at an NFW profile holds at such small radii 
[sL I25 n30| . Moreover, there are many other phenomena 
that could significantly affect the density structure of ha- 
los on these scales [23 . Effects that rely on the presence 
of large baryonic components in the satellite halos are 
probably negligible for the vast majority of subhalos that 
are otherwise dark. Other effects, such as that of tides 
and heating due to rapid encounters with the disk of the 
host and/or other subhalos are likely to be more impor- 
tant In these cases, the net result is generally 
to make subhalos more diffuse in their centers than our 
model predicts, so we feel that our prescription gives the 
largest justifiable overestimate for the gamma-ray lumi- 
nosity of subhalos. Nevertheless, due to the uncertainties 
in the density structure of subhalos at very small radii, in 
the following Sections we investigate the effect of adopt- 
ing core radii given by Tc = /3rc,o, where rc,o is the core 
radius given by the solution to Eq. (|2Jl. We allow f3 to 
vary between 10~^ and 10®. 

Figurenshows the relevant substructure properties de- 
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rived from 100 realizations of MW-like host halos. In the 
left panel of Figure ^ we show the relationship between 
the quantity pi [as we discuss below, this combination 
is roughly proportional to the gamma-ray luminosity of 
the subhalo, see Eq. Hll() ] and subhalo mass for the sur- 
viving subhalos in our model. In the right panel of Figure 
n we fix the mass range to 10^ < M/M© < 5 x 10^ and 
plot the product pi as a function of the distance of the 
subhalo from the sun, d. In both panels, we also show 
the values of pi that we computed using the method 
of Ref. [l^l- Notice that our modeling, which is cali- 
brated against N-body work, suggests that subhalos are 
typically further away from the sun and less luminous 
than was assumed by Ref. . The subhalo population 
that we present here is a more faithful representation of 
the predictions of the CDM paradigm of structure for- 
mation. Yet, we still expect that all of the simplifying 
assumptions that we have made go in the direction of an 
over-estimate of the gamma-ray signal from neutralino 
annihilations. 



III. THE GAMMA-RAY SIGNAL FROM 
ANNIHILATIONS IN SUBSTRUCTURE 

In this Section, we outline our method for calcu- 
lating the emitted flux of gamma-rays from neutralino 
annihilations in MW subhalos, and compare this with 
the gamma-ray background to determine a detectability 
threshold. 



A. Gamma-rays from a subhalo 

The number of photons with energy greater than a 
threshold energy Eth, arising from neutralino annihila- 
tions in a subhalo can be written as 



The energy-dependent effective area of the detector is 
Acs{E) and tcxp is the exposure time. As discussed 
above, for r < Tc, the density is constant, given by 
p{r) — p{rc) = Pc = const., while for r > rc, p{r) is 
given by Eq. 

The neutralino mass is M^, the thermally- averaged 
cross section times the relative velocity of the annihi- 
lating neutralinos into a final state "i" is (ajwDi, and 
dN^^i/dE is the number of photons in the energy inter- 
val E to E + dE that result from neutralino annihilations 
into final state "i." The summation runs over the differ- 
ent final states that lead to gamma-rays. 

Neutralino annihilations may yield photons in three 
ways: (1) by the direct annihilation into a two-photon 
final state (77); (2) by the direct annihilation into a pho- 
ton and a boson {jZ'^); and (3) through annihila- 
tion into an intermediate state that subsequently decays 
and/or hadronizes, yielding photons (h). The first two 
processes are easy to deal with because the result is a 
mono-energetic line emission at energies of E'^^'J = My- 

for the 77 final state and E'^^f = M^[l - {Mzo/2Myf] 
for the 7^° final state. In these instances, 
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where 

.^halo(£') 



1 r^'- d^Snsy{E) 



with n^~i = 2 for the 77 final state, and n^^o = 1 
for the jZ'^ final state. The theoretical cross sections 
for these processes span many orders of magnitude and 
vary depending upon the composition of the neutralino. 
In the context of the constrained MSSM, the highest 
values these cross sections attain are ~ 10~^*cm^s^^ 
when the neutralino is primarily composed of higgsi- 
nos, while mostly gaugino neutralinos have the lowest, 
roughly ~ 10~'^^cm'^s~^. In order to maximize the prob- 
ability of detecting gamma-rays from halo substructure, 
we choose the most optimistic set of SUSY parameters, 
thus we fix the cross sections for annihilation directly into 
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■Ci,^ioiE)AesiE)texpdE, (3) photons to {a\v\)^^ = {cr\v\)yZO = 10" 
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ydy y^^^gjd(LOS) (4) 
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contains all of the information about the subhalo, and 

dS^Susy _ 



values of M^- 

The much more complicated case of decay and 
hadronization of the annihilation products leads to 
continuum emission. The resultant spectrum for the 
most important annihilation channels can be well- 
approximated by 



dE 



^ dN^^, (a\v\) 
^ dE Ml 



(5) 



dE 
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E ' 


My 




exp 





(7) 



contains all the particle physics information. In Equa- 
tions IHl)-®, V is the distance to the differential line-of- 
sight (LOS) element over which the integral is performed. 
The energy-dependent angular resolution of the detec- 
tor is (Jg{E) and d is the distance to the center of the 
subhalo. The quantity ym.ax{E) represents a projected 
length scale on the subhalo defined by the angular reso- 
lution of the detector and given by j/maxl^-) = cre{E)d. 



where (01,03) = (0.73,7.76) bT_WW and Z'^Z'^ final 
states, (ai, Q!2) = (1.0, 10.7) for 66, (ai, 02) = (1.1, 15.1) 
for ti, and {ai,a2) = (0.95,6.5) for uu. The cross sec- 
tions associated with these processes also span many or- 
ders of magnitude, from roughly a few x lO-^^cm"^ s~^ 
to as little as six orders of magnitude smaller. Again, 
we choose the most optimistic possible value and fix the 
cross sections to all of these final states to (tT|z;|)?i = 
5 X 10~^^ cm'^s"^, independent of the neutralino mass. 
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B. Background photons 

In order to detect the signal from neutralino annihila- 
tions in substructure, it must be significantly larger than 
the contaminating noise of the gamma-ray background. 
There are two contributions to this background: (1) the 
observed background due to cosmic ray hadrons and elec- 
trons which will contribute to the background for atmo- 
spheric Cerenkov telescopes (ACTs); and (2) the difi^use 
gamma-ray background supposedly coming from astro- 
physical sources. The total contribution to the gamma- 
ray count from background photons above the threshold 
energy Eth is 



dTVB 
dE 



■dE, 



where the number of observed background photons per 
unit energy interval can be written as 



^^2.67$B(i?) ' 
with 



GeV 
(9) 



<S>b{E) = 6.4 X 10" 



1.4 X 10" 



E 



GeV 

E 

GeV 



-3.3 



+ 1.8 



E 



-2.75 



GeVy 



-2.20 



(10) 



The three background contributions to $5 are as fol- 
lows. The first term corresponds to the electron back- 
ground the second term is due to hadronization of 
cosmic rays [s^, and the third term accounts for the 
diffuse gamma-ray emission measured by the Energetic 
Gamma-Ray Experiment Telescope (EGRET) All 
three of these background contributions are relevant for 
AGTs, while for space-based gamma-ray detectors the 
only relevant background is the diffuse background due 
to astrophysical sources. A fourth source of background 
photons can be the diffuse emission from neutralino an- 
nihilations in the smooth component of the MW halo it- 
self; however, the flux in this case is more than an order 
of magnitude smaller than all the aforementioned back- 
grounds and we therefore neglect it in this work. 



C. The Detectability of Substructure 

We now discuss the conditions for detecting substruc- 
ture through neutralino annihilations and the probabil- 
ity for observing such a signal. The background count 
follows Poisson statistics, so it exhibits fluctuations of 
amplitude ^ ^/Nq. The quantity called the significance 
S = Ns/ ^/Nb therefore describes the likelihood of misin- 
terpreting fluctuations in the background as the desired 
signal from neutralino annihilation. In order to minimize 
the likelihood of counting background fluctuations, one 
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(8) FIG. 2: Left: The significance S{> £th,A/x) for VERI- 
TAS as a function of threshold energy for neutralino masses 



Mx = 100 GeV (long-dashed), = 500 GeV (solid) and 
= 5TeV (dot-dashed) neutralinos normalized to the peak 
value of the — 500 GeV neutralino case. Right: The sig- 
nificance for the GLAST experiment for masses = 40 GeV 
(solid), = 100 GeV (dashed) and = 500 GeV (dot- 
dashed), normalized to the — 40 GeV case. For VERI- 
TAS, the significance is maximized at the instrumental thresh- 



old of Ki 



50 GeV for a 



500 GeV neutralino, 



while for GLAST it is maximized at an energy threshold of 
-Eth ~ 3 GeV, for a neutralino of mass — 40 GeV. 



must demand that the significance be greater than some 
predetermined requirement, S > v. We adopt the very 
liberal detection criterion — 3. Note that both the 
background and the source photon counts are directly 
proportional to the exposure time, so that S cx ^yt^xp ■ 

The next issues that need to be addressed are the opti- 
mal threshold energy for observation and the neutralino 
mass at which the probability for detection is maximized. 
These are important considerations because the signifi- 
cance is a function of threshold energy and neutralino 
mass for a fixed exposure time, subhalo distance, and 
subhalo structure, S = S{> Eth, M^). For example, at a 
fixed neutralino mass, the function diRsusy/d-E is gener- 
ally a decreasing function of energy while at fixed energy 
it is a decreasing function of the neutralino mass. For 
AGTs and space-based gamma-ray detectors, Acs{E) in- 
creases with energy while ag{E) decreases with energy. 
The background contribution <I>b (E) is a decreasing func- 
tion of energy. We would like to investigate the best- 
case-scenario for detection, so we choose to observe at 
a neutralino mass and above a threshold energy where 
the function S{> Eth, M^) is maximized. Generally, 
S{> Eth, M-^) peaks in the continuum part of the spec- 
trum. As a result, it is most fruitful to look for evidence 
of the continuum signal and then use follow-up observa- 
tions of the line-emission to confirm that the signal is, 
indeed, due to neutralino annihilations. 

For definiteness, we consider specific examples for max- 
imizing the detectability of neutralino annihilations in 
subhalos by both an AGT and a space-based gamma-ray 
detector. For the ACT, we adopt the specifications of the 
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Very Energetic Radiation Imaging Telescope Array Sys- 
tem (VERITAS) HI. In this case S, and thus the hke- 
hhood of detection, is maximized for a = 500 GeV 
neutrahno observed above the VERITAS energy thresh- 
old of Eth — 50 GeV. We adopt the specifications of the 
Gamma-ray Large Area Space Telescope (GLAST)^ |^ 
for the space-based gamma-ray detector. For this exper- 
iment, S is maximized for a neutralino with mass just 
above the present experimental bounds js^f, « 40 
GeV, for gamma-rays observed above an energy thresh- 
old of i?th = 3 GeV. Other choices of and threshold 
energy result in smaller values of the significance. 

In Figure |2 we show the detectability condition nor- 
malized to our choices for the maximum probability of 
detection. The flattening (VERITAS) and eventual drop- 
off (GLAST) in the significance at low energies is due to 
the combined effects of the decreasing effective area of 
the detectors as well as an increased background photon 
flux. The lower significance at higher neutralino masses 
results from the suppression of flux due to the M~'^ term 
in Equation |(SJl. In principle, a similar set of values of the 
threshold energy and neutralino mass can be obtained for 
any detector given its design characteristics. 



IV. RESULTS 

Our first results deal with the likelihood of a detec- 
tion by VERITAS. Figure O shows the result of our cal- 
culation for the most optimistic case of using VERI- 
TAS to detect substructure, namely, a neutralino mass 
of = 500 GeV and gamma-rays observed above the 
VERITAS threshold energy of Eth = 50 GeV. We adopt 
a detection criterion of 5 > 3 and assume a very gener- 
ous exposure time of tcxp — 250 hours. In this Figure, 
we show the cumulative number of visible subhalos in 
the whole sky as a function of an adopted lower cut-off 
in the subhalo mass M^i^. In practice, we expect that 
there should be some cut-off in CDM substructure at 
some small mass, but any such limit is unknown and it 
is expected to be well beyond the regime where N-body 
simulations can probe. Notice that the number of ob- 
servable subhalos rises only very slowly as the minimum 
mass threshold is reduced. 

This relative flatness in the number of detectable ob- 
jects as a function of Mmin can be understood from some 
simple scaling arguments. The luminosity of a subhalo 
C, can be roughly approximated as the number of annihi- 
lations per time interval in a volume of some radius j/max 
(defined through the resolution of the detector) centered 
at the center of the subhalo. Under this approximation, 
there will be a contribution to the luminosity coming 
from the core and a contribution from the rest of the 
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volume, i.e., 

L - pM + P^d'r ^ C plrl (11) 

where C — r[. — y^j^x ^ constant that depends on the 
core radius through r[. — Vc/rg + \/(\ + rc/rsY and j/max 
through y'^^^ = 1/(1 -h ymax/^s)^- Our model yields 
a relationship between luminosity and subhalo mass, 
L ~ M°-^^, that is close to what one would expect just 
from the simple scalings ps ^ c^, rg M^^^c~^, and 
Cv ~ M^^''^'^ , which is the approximate scaling of con- 
centration with mass given by the BOl model on the 
small mass scales that we consider in this calculation 
[ssf . Knowing how the luminosity of a subhalo scales 
with mass allows us to derive a scaling relationship for the 
number of visible halos above some flux threshold. The 
number of objects visible above a particular threshold is 
given by the number of objects that are contained within 
a volume of radius d ' , such that the observed flux is 
higher than the imposed threshold. In this case, the num- 
ber of objects above this threshold per logarithmic mass 
interval is dN{Ng > ^V^)/dlnM ~ (dn/dlnM) L^/^ 
Simulations, as well as our simgle, semi-analytic model, 
show that dn/dlnM - Af""** illil- Using the result 
that L ~ M"-^^, the number of visible halos per logarith- 
mic mass interval is then dA^(> vy/NB)/dlnM - M-" 02 
on small mass scales. The dependence on minimum mass 
cut is weak: making the mass cut small does not neces- 
sarily lead to a marked increase in the number of visible 
subhalos. 

We find that with a lower mass cut at Afmin — 10^ M©, 
we expect A^totai ~ 17 subhalos to be detectable at > 3. 
As we discussed above, considering lower mass cuts does 
not increase the number of detectable objects apprecia- 
bly. Even in this best-case-scenario, the 68% range only 
extends up to iVtotai — 22 (so 84% of the realizations had 
fewer detectable subhalos) and the 95% range extends 
to iVtotai — 25. This means that, on average, a detector 
like VERITAS would have to survey 1 /20 of the entire 
sky to find one subhalo; this amounts to one subhalo per 
- 170 VERITAS fields of view. Considering the smaU 
field of view of VERITAS and the fact that ACTs can, 
on average, observe only ~ 6 hours per night, such a de- 
tection must rely heavily on serendipitous discovery, even 
under the most favorable of circumstances. Note that for 
neutralino masses higher or lower than ^ 500 GeV, 
the total number of subhalos detectable at S* > 3 would 
be smaller. 

If the direct detection of substructure via particle an- 
nihilations (or lack of such a detection) is to yield any 
constraint on SUSY and/or information about struc- 
ture/galaxy formation, it is necessary to investigate the 
infiuence of the background cosmology on the expected 
number of visible subhalos. As an example, in Figure 
13 we show the results of a calculation based on a cos- 
mology with Qm = 1 ^ = 0.3, but with a nonstan- 
dard primordial power spectrum of density fluctuations. 
In this particular case, we consider a power spectrum 
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FIG. 3: The cumulative number of subhalos with mass M > 
Afmin that are detectable at a significance S > 3 by VERITAS 
on the entire sky as a function of Mmin. As discussed in the 
text, the neutralino mass is = 500 GeV observed above 
the VERITAS threshold energy of Eth = 50 GeV, yielding 
the highest possible detection efficiency. The exposure time 
is set to a generous iexp = 250 hours. The solid line represents 
the mean over all model realizations, down to the minimum 
mass threshold of 10^ M©, in the standard ACDM cosmolog- 
ical model with scale-invariant primordial power spectrum. 
The error bars demarcate the 68 percentile range (symmetric 
about the median) determined from the model realizations. 
The dashed line represents the same region for a standard 
ACDM model with a "running" power law spectrum of pri- 
mordial fluctuations with dn/dlnk = —0.03 as suggested by 
the WMAP team The down arrows indicate that more 
than 16% of the realizations had zero detectable subhalos in 
the corresponding mass bin. 



with a running power law index dn/dlnk = —0.03, with 
n{k = 0.05 Mpc"^) = 0.93, and as = 0.84, as suggested 
by the recent analysis of the WMAP data l] . Note that 
the statistical significance of this result is weakened if 
additional uncertainties in the mean flux decrement in 
the Ly-a forest are considered jS^; nevertheless it is of 
interest to investigate the consequences of such a power 
spectrum. In this alternative model the mean number 
of visible subhalos is reduced by a factor of ~ 50. This 
effect can be understood by examining the influence of 
reduced small-scale power on the average properties of 
CDM substructure 0, 113 ■ In models with reduced 
small-scale power, low-mass objects form later, and as a 
result have lower concentrations than objects of the same 
mass in the standard model. The reduced concentrations 
of typical subhalos and the fact that the subhalos are 
consequently more susceptible to tidal disruption, and 
therefore typically exhibit a larger core in their radial 
distribution within the host, results in a large reduction 
in the number of potentially detectable subhalos relative 
to the predictions of the standard model. 



FIG. 4: The cumulative number of subhalos on the entire sky, 
with mass M > Mmin, that are detectable at a significance 
S > 3 by GLAST as a function of Mmin. We show results for 
two choices of neutrahno mass, = 40 GeV (solid lines) 
and — 100 GeV (dashed lines), observed above an en- 
ergy threshold of Eth = 3 GeV. The error bars show the 68 
percentile range. As in Figure |3 down arrows indicate that 
more than 16% of the realizations had zero visible halos in 
the corresponding mass bin. 



We now turn our attention to the detectability of sub- 
structure with GLAST. In Figure 0] we exhibit the num- 
ber of subhalos that are detectable at 5 > 3 after a year 
long exposure with GLAST. Here the most optimistic 
number of detectable subhalos is A^totai ^ 19 at the upper 
limit of the 68% range. This amounts to ^ 2 objects per 
field of view, on average; however, it is necessary to exer- 
cise caution. Consider the energy scales involved in this 
calculation. The optimum neutralino mass for a GLAST 
detection is at the lower limit of current experimental 
bounds, ~ 40 GeV [13 . The factor dSRsusy/di; in 
Eq. Q is proportional to M~'^ and because of its small 
effective area, GLAST cannot compensate for this rapidly 
decreasing function of M^. Consequently, the number of 
halos that are observable drops dramatically as is 
increased. This is shown by the bottom set of (dashed) 
lines in Fig. 0| 

Due to the uncertainty in the density structure of the 
inner regions of halos, it is of interest to investigate the 
sensitivity of our results to the choice of core radius. In 
Figure [S] we show how our results on the number of de- 
tectable subhalos vary with core radius. We parameter- 
ize the variation in core radius by (3 = Tc/rco, where 
rcfi is the core radius assigned according to our standard 
procedure of solving Eq. (j^J , and Vc is the new core ra- 
dius defined as a multiple of r^o- Clearly, the particular 
choice of the core radius affects our results only weakly as 
a shift in core radius of nine orders of magnitude changes 
the number of detectable subhalos by less than 40%. This 
robustness is easily understood. As we discussed previ- 
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FIG. 5: The cumulative number of detectable subhalos with 
mass M > Mmin as a function of Mmin for different values 
of the quantity P — fc/fcfi (see text). The neutralino mass, 
exposure time, threshold energy and instrument specifications 
are as in Fig. |3 Note that over a change of nine orders of 
magnitude in /3, the total number of visible halos is nearly 
independent of the core size. 



ously, the luminosity of a subhalo can be approximated 
by L Ps'^s (''c~2/max) [^^6 Eq. lfTT|l ]. For a fixed subhalo 
mass at a fixed distance (i.e., y^n^x — const.), « 1 as 
long as Tc << Ts and so the expected luminosity remains 
approximately constant as the core radius is changed. 
This behavior is valid up to approximately /3 ~ 10^. At 
this point, starts to deviate from 1, the core size comes 
nearer the typical range of ymax values, and the constant 
C decreases. Given the spatial distribution of the sub- 
halos in our model, this behavior becomes noticeable at 
roughly Tc/ts ^ . These scales are beyond the cur- 
rent resolution limits of the inner parts of subhalos in 
N-body simulations. In this case, our results suggest 
that the presence of smaller cores does not necessarily 
mean that the number of potentially detectable subhalos 
increases appreciably. 



V. DISCUSSION AND CONCLUSIONS 

We considered the detection of otherwise dark CDM 
substructure from the gamma-ray products of neutralino 
annihilations in the center of subhalos. In order to 
demonstrate several important points, we chose the most 
optimistic SUSY parameters so that we maximized the 
probability of detection. We also employed a realistic, 
yet still optimistic from the standpoint of predicting ob- 
servable signals from substructure, model for the popu- 
lation of subhalos in the MW. Our analysis allowed us to 
estimate the likely range in the number of detectable sub- 
halos, given the statistical fluctuations in subhalo popu- 



FIG. 6: The number of subhalos that are detectable by VER- 
ITAS at S > 3 as a function of the minimum mass cut-off 
Mmin, after adopting a Moore profile rather than a NFW 
profile to describe the substructure (short-dashed) or allow- 
ing a positively tilted primordial power spectrum with n = 1.1 
and COBE normalized to erg — 1.2 (long-dashed). The solid 
line represents the prediction of our fiducial "standard" model 
of NFW halos and a scale-invariant power spectrum with 
CT8 — 0.95. The error bars are as in Figure |3 

lations from host to host. The simplifying assumptions 
of our model should only lead to overestimates of the 
core densities of subhalos and underestimates of the typ- 
ical distance to the nearest subhalo and, if anything, we 
expect that the actual number of detectable subhalos 
should be smaller than the results that we presented here. 

We found it unlikely that an ACT similar to VER- 
ITAS will be able to detect neutralino annihilations in 
dark subhalos due to the small number of potentially de- 
tectable subhalos and the relatively small field of view of 
such a detector. On the other hand, we found that if the 
mass of the neutralino is near the low end of the viable 
range, 75 GeV, the GLAST detector, with its much 

larger field of view (~ 2 steradians) will have a signifi- 
cantly higher probability of detecting such objects. In our 
best-case-scenario of a neutralino of mass = 40 GeV 
and an observing threshold of Eth — 3 GeV, we found, on 
average, two subhalos visible at S* > 3 per GLAST field 
of view. Our results indicate that the best strategy is to 
locate gamma-ray signals with GLAST and then point 
detectors like VERITAS at the location of these signals in 
order to pinpoint them, observe the shape of the gamma- 
ray spectrum, and search the energy spectrum for the 
line-emission at = M^, which is a "smoking gun" of 
neutralino annihilation in substructure. As an example, 
we found that for a neutralino of mass ~ 75 GeV, 
there will be ~ 1 detectable subhalo per GLAST field of 
view, on average, assuming the maximum cross section 
for annihilation into photons, and that a directed VER- 
ITAS exposure should confirm the line-emission feature 
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after an exposure time of t^xp ^ 450 hr. If the neutralino 
mass were to be in the range, 100^ M^/GeV^ 500, we 
found that a detection would require an instrument with 
a large effective area, such as VERITAS; however, such 
a detection would have to rely on serendipity because 
of the comparably small field of view of such detectors. 
Unfortunately, for M^> 500 GeV we found that it is un- 
likely that neutralino annihilations in dark subhalos will 
be detectable. 

As we mentioned earlier, our results are rather sensitive 
to the assumed form of the subhalo density profile. To 
illustrate the significance of this, we have repeated our 
calculations assuming the density profile of Moore et al. 
P, p{r) = Po{r/rM)-^/Hl + rlruY^^^. with p(r) cx 
r~^l'^ at small radii. Our results are presented in Figure 
El Adopting the steeper Moore et al. profile increases 
the number of detectable subhalos by approximately an 
order of magnitude. If subsequent studies reveal that 
CDM predicts density profiles that vary as pir) oc r~'', 
with 1 < 77 < 1.5 intermediate between the Moore et 
al. and NFW profiles, then the number of detectable 
subhalos should lie between the solid and short-dashed 
lines in Fig. 

Suppose that gamma-rays are not observed from halo 
substructure. Our results demonstrate that the lack of 
such a detection doesn't lead to a bonanza of constraints 
on SUSY in general, or the constrained MSSM in partic- 
ular, because the number of subhalos that are detectable 
via gamma-rays depends sensitively on the properties of 
the substructure. Even after choosing the optimal pa- 
rameters for detection (i.e., the largest couplings to pho- 
tons) the likelihood of a detection is small for most of the 
viable mass range of the neutralino. Additionally, as we 
have shown, the likelihood depends upon the density pro- 
files in the innermost regions of dark matter halos. The 
work of Power et al. |30,J shows no evidence that density 
profiles approach the NFW inner power law of p oc r^^ . 
Rather, they find that the power law may become ever 
shallower with decreasing radius. Moreover, there are 
additional uncertainties that are not associated with our 
lack of knowledge of density profiles and subhalo pop- 
ulations. The predictions of the number of potentially 
detectable subhalos are strongly dependent upon poorly- 
constrained cosmological parameters. We showed in Fig. 
I^lthat adopting the best-fitting power spectrum from the 
WMAP group reduces the probability of detection by a 
factor of ^ 50 relative to the same cosmological model 
with a standard, scale-invariant primordial power spec- 
trum. 

Of course, a detection of gamma-rays from substruc- 
ture would yield a plethora of information. First and 
foremost, it would be strong evidence for neutralino (or 
some other WIMP which annihilates into a photon fi- 
nal state) dark matter, it would indicate the presence 
of numerous, otherwise dark subhalos surrounding the 
MW and disfavor models with reduced small-scale power 
[lol[llj |. and it would indicate that such subhalos achieve 
high densities in their central regions. If future detec- 



tions yield surface brightness profiles for substructure, 
then the dark matter distribution of the subhalos could 
be inferred which would be an important addition to the 
studies of mass density profiles [23, 113, 113, 113 • If ^ large 
number of subhalos are found then that would give us 
information on the survival of the very inner regions of 
accreted subhalos and provide an insight into the accre- 
tion history of the MW halo. However, we must exercise 
caution: ij substructure in the Milky Way halo is found 
via annihilation into high-energy photons, it will still be 
difficult to "measure" supersymmetric parameters. First, 
as we have already mentioned, the expected number of 
detectable subhalos is sensitive to cosmological parame- 
ters. In Fig. El we show the number of subhalos that are 
detectable by the VERITAS instrument in the case where 
the power spectrum has a "blue tilt," with n = 1.1, and is 
normalized to the Cosmic Background Explorer measure- 
ments of cosmic microwave background anisotropy |60j| . 
yielding erg — 1.2. In this case the SUSY parameters are 
fixed, yet the number of detectable halos is boosted by 
a factor of ~ 3 due to the increase in small-scale power. 
Although spectra with such large values of erg seem to 
be disfavored (e.g., H^), uncertainties like these must 
be marginalized over in analyses of the SUSY parame- 
ter space and calculations similar to the one presented 
here may be helpful in this regard. Second, the observed 
flux of gamma-rays from a particular subhalo depends 
upon the distance to that particular subhalo. This is a 
significant problem, for there is no obvious way to deter- 
mine reliably the distance to an otherwise dark subhalo. 
In our approach, we have attempted to marginalize over 
such uncertainties by calculating "likely" realizations of 
substructure in the MW and the most likely number of 
observable subhalos, assuming a maximum annihilation 
cross section to photons. Therefore, if a large number of 
subhalos are observed, it is likely that the neutralino an- 
nihilation cross section is large. Such a large cross section 
would imply a neutralino with a large higgsino compo- 
nent, a region of the supersymmetric parameter space 
which will be more easily explored in forthcoming exper- 
iments. If, on the other hand, SUSY is discovered in the 
upcoming runs of the Large Hadron Collider^ (LHC), 
and the mass and couplings of the neutralino are sub- 
sequently measured, this information would provide the 
framework within which a calculation along the lines of 
the one presented in this manuscript could be used to 
study the properties and distribution of CDM substruc- 
ture. 
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